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ABSTRACT 

We use the Direct Simulation Monte Carlo (DSMC) method combined with an n- 
body code to study the dynamics of the interaction between a gas-rich spiral galaxy 
and intracluster or intragroup medium, often known as the ram pressure scenario. 
The advantage of this gas kinetic approach over traditional hydrodynamics is explicit 
treatment of the interface between the hot and cold, dense and rarefied media typical 
of astrophysical fiows and the explicit conservation of energy and momentum and the 
interface. This approach yields some new physical insight. Owing to the shock and 
backward wave that forms at the point ICM-ISM contact, ICM gas is compressed, 
heated and slowed. The shock morphology is Mach- disk- like. In the outer galaxy, the 
hot turbulent post-shock gas fiows around the galaxy disk, while heating and ablating 
the initially cool disk gas. The outer gas and angular momentum are lost to the fiow. 
In the inner galaxy, the hot gas pressurizes the neutral ISM gas causing a strong 
two-phase instability. As a result, the momentum of the wind is no longer impulsively 
communicated to the cold gas as assumed in the Gunn-Gott (1972) formula, but 
oozes through the porous disk, transferring its linear momentum to the disk en masse. 
The escaping gas mixture has a net positive angular momentum and forms a slowly 
rotating sheath. The shear fiow caused by the post-shock ICM fiowing through the 
porous multiphase ISM creates a strong Kelvin-Helmholtz instability in the disk that 
results in Cartwheel-like ring and spoke morphology. 
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1 INTRODUCTION 

The interaction between the intracluster medium (ICM) or 
the intragroup medium (IGM) and the interstellar medium 
" (ISM) in a gas-rich spiral is often identified as the source 
' of star bursts and gas starvation followed by morphological 
transformation. Recent observations of this process in the 
Virgo cluster have spawned a number of detailed hydrody- 
namic simulations that investigate the conditions required 
for these scenarios. This paper applies the hyb rid Direct 
Simul ation Monte Carlo (DSMC) method from [Weinber g 
(|2013l , hereafter Paper 1) to the collision of two distinct 
gas components with different densities, temperatures and 
bulk velocities. Such an interaction is typically referred to 
as ram pressure. Many if not most ram pressure interactions 
are in the transitional regime which is not quite hydrody- 
namic. The subsequent evolution of the system depends on 
the dynamics at the collision interface and thereby provides 
a good test of the kinetic approach. Thus, the advantages 
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of DSMC all come to the fore when considering the ram- 
pressure mechanism and make ram pressure both a good 
problem for code comparison. 

Although DSMC was developed for rarefied and transi- 
tional gas fiows, the method has advantages for near con- 
tinuum (i.e. low Knudsen number) fiows: 1) there is no 
Courant-Friedrichs-Lewy (CFL) condition for DSMC; that 
is the method is unconditionally stable; 2) the method is 
100% shock capturing; and 3) the DSMC approach naturally 
resolves the gas-phase interfaces that plague many hydrody- 
namic simulations when applied to scenarios with ISM cool- 
ing, e.g. giving rise to the over-cooling problem. The penalty 
for these advantages is much lower computational efiiciency 
than the same computation using traditional grid hydro- 
dynamics. This approach requires many more gas particles 
in the continuum limit than mesh cells or SPH particles, 
although it has the advantage of consistently transitioning 
from dense to dilute gases, fully capturing shocks, and re- 
solving phase boundaries. By design, DSMC is correct even 
as the mean free path becomes large. Such regimes are com- 
mon in intracluster gas dynamics, hot bubbles and conduc- 
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tion at interfaces, and in multiphase ISM on small scales 
which require accurate treatment of interfaces. Moreover, 
this method maintains the underlying physics of the multi- 
phase interfaces. 

As the work proceeded, it became clear that the stan- 
dard explanation of ram pressure did not explain results 
of the numerical experiments, either those described here 
or published elsewhere. To illustrate these differences, we 
will begin with the standard explanation for astrophysical 
ram pressure in ^ followed by an exploratory sequence of 
simulations of increasing generality in JSl These results will 
be compared with published observational and theoretical 
findings in 311 In brief, many of the features found here 
have been identified and reported by others (as expected 
for a code test) but the simulations here reveal a com- 
plexity of multiphase instabilities and angular momentum 
transport in the cold gas disk that has not been described 
previously. The overall explanation most commonly offered 
for this interaction is th e ram-pressure paradigm envisioned 
bv lGunn fc GottI (|l972l , hereafter Gunn-Gott). However, we 
will argue that the dominant processes driving the evolu- 
tion of a galaxy with a wind is not Gunn-Gott ram pressure 
but the combination of shock heating, boundary- layer abla- 
tion, Kelvin-Helmholtz instability in a multiphase medium, 
and angular momentum transfer. Because of this complex- 
ity, any specific observed interaction is likely to depend on 
the details of the galaxy, its gravitational environment and 
the density and speed of the wind. Owing to the expense 
of the DSMC method, we can not afford an exhaustive sur- 
vey, but choose a single galaxy model and a high-density 
wind of modest to low velocity. In particular, we explore 
and describe the results for a one-dimensional slab, a three- 
dimensional disk + halo model with face-on, oblique, and 
edge-on winds, and include a case with and without a hot, 
coronal halo. We summarize in JSj 



2 RAM PRESSURE 

2.1 The astronomical setting 

Galaxy clusters and groups are often permeated by hot, ten- 
uous. X-ray emitting gas known as the intracluster medium 
(ICM). Throughout this paper, we will refer to both the in- 
tracluster medium and the intragroup medium (IGM) as the 
ICM. The pressure exerted by flows in this medium on a cold 
gas disk bound to the stellar disk is often described as ram 
pressure. In physics and engineering, ram pressure denotes 
the back pressure applied to a solid body moving through a 
fluid: P — pv where P is the pressure, p is the mass density, 
and V is the relative speed of the body. Astronomers use ram 
pressure to denote the the pressure applied by one fluid on 
another when the two fluids have distinct phases and sig- 
nificant relative bulk velocity. Let us refer to a hot tenuous 
gas that has significant bulk velocity with respect to the 
galaxy disk as a wind. Ram pressure stripping is said to oc- 
cur if this wind is able to impart sufficient momentum to the 
cold disk that it overcome the gravitational attraction of the 
galaxy. Although we will retain the traditional terminology, 
this paper will demonstrate that ram pressure is a physical 
misnomer in the astrophysical context; momentum transfer 
does occur, but it does not act like a macroscopic blunt body 



moving in a hydrodynamic flow. Rather, the details of the 
interaction of the two gas components with differing densi- 
ties and temperatures is sufficiently dynamic that the ram 
pressure paradigm misses may of the essential features. 

Regardless of the dynamical details, observational evi- 
dence for the interaction clearly exists. For simplicity, let us 
follow the standard practice of referring to the flow in hot 
medium relative to the frame of the cold galaxy as a wind 
and the resulting momentum transfer as ram pressure. Evi- 
dence for ram pressure can be found in many galaxy clusters. 
For example, three morphological features in NGC 4404J 
suggest strong gas interactions. First, the bowed morphology 
of the gas suggests a balance between the ram pressure and 
the lower gravitational force at larger disk radii. Secondly, 
the lower gas content in the extended blue stellar disk sug- 
gests a recent epoch of star formation, presumably induced 
by the interaction event. Finally, lanes of reddening can be 
seen trailing the inferred motion of galaxy in the cluster. 
The expected result of ram pressure stripping is a galaxy 
which contains very little cold gas. This effectively halts 
star formation in the galaxy, supporting the belief that ram 
pressure stripping could be one of the pr ocesses responsible 
for th e morphology-density re lation (see van der Wei et al.l 
2010, and references therein). ICrowl et al.l (|2007|) describe 
the observational case for ram-pressure stripping of NGC 
4402 in detail. 

Another interesting candidate is NGC 4524j. Line-of- 
sight velocities suggest that this galaxy is currently falling 
towards the center of the Virgo cluster (toward the bottom 
of the HST image^). The wind direction is less certain but 
presumably from below. The bowed and truncated disk, and 
the concentration of dust and gas to one side of the galaxy 
are all indicators that ram pressure stripping is forcing gas 
out of the galaxy. H-alpha images of this galaxy suggest 
that there is no gas in the outer disk, but there are gas 
fllaments emerging from the plane of the disk, presumably 
caused by the interaction with the intracluster medium. The 
plume geometry of these features in NGC 4522 is consistent 
with predictions by the simulations described in ^3.41 that 
post-shock gas will excite and then ooze through multiphase 
instabilities in the cold disk. 



2.2 The Gunn-Gott ram-pressure criterion 

We review the physical arguments for the classic Gunn-Gott 
ram pressure formula before moving on to the numerical 
experiments using DSMC in the next section. Let us begin 
by considering a parcel of wind gas moving with respect the 
cold disk. The momentum per unit volume in the wind is 



wind momentum 



■ pICMVwind- 



(1) 



The wind parcel crosses the stellar and gaseous disk in the 
time interval approximately given by h/vwind- Therefore, 
the momentum delivered by the stream of ICM parcels to 
the cold gas disk is 

wind force = picMV^ind x Vwind/h . (2) 

l/duration 



Hubble image: 
Hubble image: 



http://www.spacetelescope.org/images/heic0911c 



http://www.spacetelescope.org/images/heic0911b 
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The horizontally bracketed term in the equation above is the 
arrival frequency of ICM gas parcels, ignoring any dissipa- 
tion. Finally, the condition that the volume force applied by 
the wind exceeds the gravitational restoring force binding 
the cold gas to the stellar disk is 



2 /^ ^ ^^ 

PlCMVy,ind/ll > Pgas-^. 



(3) 



Rearranging, we get a condition on the ram pressure re- 
quired to strip the disk gas: 



2 . ^ 

■ PiCMVwind > ^gas 



dz 



(4) 



3 EXPERIMENTS 

We will attempt to verify this physical picture with a short 
series of numerical experiments of increasing complexity. 
The numerical methods described in Paper 1 are briefly re- 
viewed in the first section below and the initial conditions 
for the experiments are presented in ^3.21 These are followed 
by in-depth descriptions of the one-dimensional f ^3.3|) and 
three-dimensional ( ^3.4|) hot wind-cold disk simulations. 



3.1 Methodology 

DSMC solves the collisional Boltzmann equation when the 
collision term is dominated by two-body collisions for one 
or more species with or without internal degrees of free- 
dom. Unlike molecular dynamics, the solution is obtained 
by sequentially computing the collisionless or ballistic flow 
followed by simulating the collisions implied by the Boltz- 
mann collisional term. Essentially, the latter step is a Monte 
Carlo evaluation of the collision integral, and therefore a sin- 
gle simulation is one realization of an ensemble of possible 
solutions. To help prevent the near-continuum limit from be- 
ing a computational bottleneck. Paper 1 describes a hybrid 
DSMC approach that recovers the Navier-Stokes equation 
limit for mean free paths that are very small compared to 
any flow scale. This is achieved by limiting the number of 
collisions when the mean-free path becomes very small com- 
pared to any other scale of interest. This algorithm was par- 
allelized and merged with a collisionless N-body solver called 
EXP. Paper 1 demonstrates that this code correctly repro- 
duces the standard shock-tube tests and develops Kelvin- 
Helmholtz instabilities following the analytic dispersion re- 
lation. 

In the simulations below, the stellar disk and dark mat- 
ter halo are self- gravitating and mutually interacting. The 
gas feels the gravitational field generated by the stellar and 
dark matter components but do not exert gravitational force 
on these components or on itself. Therefore, in these runs, 
the gas cannot be gravitationally unstable. This is not a fun- 
damental limitation of the combined EXP-DSMC code. For 
example, EXP also includes a PIC-Fourier Poisson solver 
that could be used to investigate the small-scale gravita- 
tional response of the gas to the momentum transfer from 
the wind. Rather, because the gravitational instabilities in 
the gas are likely to be on a much smaller scale than those 
of the stellar and dark matter components and because the 
mass in diffuse gas is very small compared to the mass in the 



stellar disk or halo, we choose to focus on the large-scale fea- 
tures of the flow. A simulation with a coronal halo is briefly 
described in ^3.61 the heat generated by the wind interaction 
is sufficient to prevent the coronal gas from cooling. 

3.2 Initial conditions 

The first set of experiments (described in detail in ^3.3p 
explores a one-dimensional wind interaction with a slab of 
cold gas confined by stars and dark matter. A planar stel- 
lar slab is realized from the exact isothermal solution to the 
equations of hydrostatic equilibrium. In the second set of 
experiments f ^3.4D , the stellar disk and h a lo are initialized 
as described in IHollev- Bockelmann et al. (2005). In brief, 
the steps in realizing the phase-space distribution for the 
second set of simulations are as follows: 1) the dark halo 
distribution function is realized b y the Eddington inversion 
method fe.g. lBinnev fc Tremaine 1987), modified to include 
the monopole component of the stellar disk. The stellar 
disk is realized using Jeans' equations with the potential of 
the disk and halo particles as computed by the EXP Pois- 
son solver. The profiles are Milky- Way like: the dark halo 
is an NEW (Navarro et al. 1996) profile with c = 15 and 
Rvir = 300 kpc and an exponential stellar disk with scale 
length of 3 kpc and scale height of 350 pc. The stellar to dark 
matter ratio is 1/15. In both cases, the cold gas disk is ini- 
tially in local hydrostatic equilibrium with a temperature of 
3000 K based on numerical solutions of the one-dimensional 
Euler equations. For the Milky- Way like model, the gas disk 
has the same scale length as the stars and is rotationally 
supported. We do not include a separate bulge or spheroidal 
component in this model. 

As described in Paper 1, these tests use an LTE cool- 
ing approximation rather than the full cross-section based 
Monte Carlo simulation of the Boltzmann collision term to 
facilita te most of t he published work on the subject. This 
model (Black 1981) includes atomic processes only so that 
temperatures lower than 3000 K will be of little energetic 
consequence. We will describe gas whose temperature is at 
the minimum of the cooling curve 'cold' gas. Use of the full 
DSMC method allows gas of any chemical composition and 
excitation state, and this will be studied in a later contribu- 
tion (see 35] for additional details). The hot halo used in ^3.61 
is realized for simulation solving the hydrostatic equilibrium 
equation for the monopole term of the total gravitational po- 
tential fixing the temperature to T = 10^ K and the density 
at the solar circle to n = 10 ~^ atom/cc. The wind is real- 
ized as a uniform column with T = 10^ K and n = 10 ~^ 
atom/cc, initially moving toward the galaxy with a uniform 
bulk velocity of 2 x Vvir = 240 km/s. This modest wind ve- 
locity is low compared to that expected in a large cluster 
investigated by some authors (see SJ- However, simulating 
a high- velocity wind requires a large particle column and is 
computationally expensive, so we compromise with a wind 
more typical of a group interaction. We explore three wind 
directions relative to the rotation axis of the disk: 0, 45 and 
90 degrees. 

The code is heavily instrumented, accumulating the 
mean Knudsen number, comparing the mean ballistic trajec- 
tory to the cell size, collisions per simulation particle, frac- 
tion of particles and cells in the EPSM limit for every cell at 
every time step. The mean velocity field, (v), is character- 
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ized by gridding the gas particles and computing streamlines 
by numerically solving the first-order differential equations 
for the flow: x = (v). Other mean quantities of the flow 
such as energy and angular momentum fluxes are readily 
computed from the gas particle distribution. These ensemble 
diagnostics will be used to illustrate the field quantities in 
the following sections. However, unlike hydrodynamic simu- 
lations, the particles in a DSMC simulation are not tracers 
of the underlying fields but carry the intrinsic quantities of 
the microscopic physical state of the gas. The computation 
of ensemble quantities, then, artificially degrades the true 
resolution of the simulation. In fact, the gas dynamics in 
the simulation takes place on a scale much smaller than the 
ensemble grid scale, often by an order of magnitude. For 
scenarios with unsteady flow, an accurate spatial ensem- 
ble estimate requires averaging over multiple simulations. 
This applies to hydrodynamic simulations as well, but this 
is rarely done in astrophysical simulations since a precise 
estimate at a point in space has little predictive value. 

Note that we will use two sets of units in this paper: 1) 
system units that assume that the mass and length of the 
entire bound system along with the gravitational coupling 
constant is unity; and 2) Milky Way units that scale the 
model to our Galaxy using physical units-kiloparsecs, years, 
and km/s. When quantities are quoted without units, they 
refer to system units and when quantities are physical units, 
they refer to Milky Way units. 

3.3 Slab model: simple experiment 

The first experiment is that of a one-dimensional slab of 
stars whose gravity confines a cool gas layer. This was de- 
signed to explicit test the Gunn-Gott criterion in its natural 
domain of validity: a wind impinging on a cold gas layer 
that is gravitationally confined to a plane-parallel slab. The 
stars are distributed with the classic isothermal distribu- 
tion p(x,y,z) oc sech^(z//i) with scale height h. The gas is 
initialized in isothermal equilibrium with the confining slab 
with T = 3000 K and mid plane density of n = 1 atom/cc. 
The boundary conditions are periodic in the plane of the 
slab (x and y directions) and vacuum in the vertical (z) 
direction. The simulation begins with a hot homogeneous 
wind (T = 10^ K, n = 0.001 atoms/cc) filling the the region 
z < —lOh with a mean vertical velocity chosen to exceed 
the Gunn-Gott criterion by a factor of three. As described 
above, all gas components may cool by atomic and plasma 
cooling mechanism following the standard LTE cooling com- 



putation. See Figure 1(a) for a schematic description. 



Much to my surprise, the cold gas does not leave disk! 
Rather the evolutionary scenario is as follows. The gas 
shocks as it approaches the disk. The shock heats the gas, 
and the upstream compression wave propagates back into 
the original wind. The velocity of the turbulent, denser cool- 
ing layer of the post-shock hot gas is lower than v^ind and 
the pressure begins to displace the cold layer from mid plane. 
The momentum from the wind is transferred to the bulk 
ISM through a cascade of primary and secondary inelastic 
scattering events. See Figure |l(b)| for the location of these 
regimes. Furthermore, the hot gas begins to pressurize the 
cold gas, exciting two-phase instabilities that give rise to 
temperature and density 'fingers' as the 10^ K gas displaces 
the volume originally occupied by the cold gas. We will see in 
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(b) the wind-disk interaction 

Figure 1. Description of the initial conditions | (a) | and the evolv- 
ing hot-wind-cold-disk interaction |(b)| in the frame comoving with 
the disk. The cold and hot gas is indicated by blue and red shad- 
ing respectively. In addition, the post-shock gas is shaded yellow 
and an enhanced cooling layer at the hot-cold interface is shaded 
purple. 



]that this production of multiphase features is enhanced 
by rotational shear. The total momentum is conserved but 
the slowed, compressed flow downstream of the shock im- 
plies that the wind does not transverse the disk on a ballistic 
crossing time. Rather, the slowed flow of hot gas increases 
the effective interaction time and allows the stellar disk to 
respond dynamically to the displaced cold disk. Then, the 
hot wind's momentum is transferred to the stellar disk grav- 
itationally. So in the end, the cold gas layer is not lost but 
rather the entire disk, gas and stars together, is accelerated. 
This simulation was repeated with gas cooling. With- 
out cooling, the energy deposited by the wind is sufficient to 
heat the gas layer to 10^ K in 200 Myr. The cold gas then 
escapes the gravitational pull of the disk. Although this re- 
sult agrees with the Gunn-Gott prediction, the mechanism 
is different: evaporation rather than impulsive momentum 



© 2013 RAS, MNRAS QQQ.nHTGl 



Ram Pressure Dynamics 5 



transfer. Cooling is an essential ingredient that prevents the 
cold gas disk from evaporating. 

In summary, the scenario described in ^2.21 failed to de- 
scribe the flow because the momentum transfer is not in the 
impulse limit as derivation of equation Q assumes. Owing 
to the shock, the post shock gas acts as a slowing moving 
piston. This gives time for the stellar disk to respond grav- 
itationally to the displacement of the cold disk from the 
stellar mid plane. As a result, the initially cold gas layer is 
decelerated by the stellar disk causing the momentum from 
the wind to be transferred to the disk, and in the long run 
the entire system, although the disk-halo interaction occurs 
on a much larger timescale than in this simulation. 

3.4 Galaxy model: harder experiment 

The first experiment was intentionally idealized and de- 
signed to check the Gutt-Gott mechanism in the simplest 
conceivable scenario. Intuitively, the inhomogeneity in the 
galaxy density profile will change both the downstream flow 
patterns and possibly entrain cold gas. How will this change 
the simple scenario above from ^3.3P 

To test this, we use the more realistic equilibrium galaxy 
model with the same wind column as in the previous ex- 
periment. The wind is only distributed in the vicinity of 
extended disk to best use computational resources as de- 
scribed in ^3.21 To prevent the adiabatic expansion of the 
wind column perpendicular to the flow, the gas is subjected 
to periodic boundary conditions in the two perpendicular 
directions. The duration of the simulation then is limited 
by either (i) the complete passage of the finite wind column 
past the galaxy or (ii) the wrapping of the flow transverse 
to the wind direction, which ever happens first. The ver- 
tical extent of the wind column is chosen, empirically, so 
that these two conditions obtain at approximately the same 
time. For simulations here, this is approximately 500 Myr 
after the wind impact, scaled to a galaxy of Milky- Way size. 

Vertical slices containing the geometric center of the 
disk through the gas density and temperature fields con- 
structed from the ensemble values of the DSMC particles 
are shown in Figures [2]— [3 The initial condition, just before 
the vertically moving wind reaches the disk is depicted in the 
top- left panels of Figures [2] and [S] The disk plane is initially 
at 2; = 0. As in the simple slab model, the wind gas shocks 
as it reaches the galaxy. Initially, the momentum of the wind 
atoms is transferred directly to the cold gas. However, after 
the shock forms, the momentum from the incoming wind is 
transferred by the shock to the slower, compressed down- 
stream flow. The heated post-shock gas slows and piles up, 
increasing its density and driving a compression wave up- 
stream into the wind. After approximately 100 Myr (the 
lower-right panels in Figs. [2] and O, a well-defined Mach 
disk has formed with a temperature of close to 10^ K, offset 
approximately 0.02 units below the disk plane (approx. 6 
kpc, scaled to the Milky Way). Meanwhile, the hot wind be- 
gins to flow around the disk, heating and ablating the outer 
initially cold disk gas. Most of the disk gas is ablated at 
approximately 3 disk scale lengths. The temperature of the 
ablated gas is approximately 8 x 10^ K. This gas cools to ap- 
proximately 4 X 10^ K as it expands and forms a bell-shaped 
sheath. 

Figures [3] and [6] describe the gas flow between 150 to 







Figure 2. Edge-on view of the combined hot- wind-cold-disk gas 
evolution. Each panel describes the density in units of hydrogen 
atoms per cubic centimeter. One length unit is approximately 300 
kpc scaled to the Milky Way galaxy. The elapsed time is inset in 
the lower right-hand side of each panel. The flow is cylindrically 
symmetric on large scales. Fig. [3] 







Figure 3. Continuation of Fig. [2] The cold, dense disk is seen 
between z = —0.01 and 0.01 (wind direction) and x = —0.04 and 
0.04 (transverse direction). The notches of lower density in the disk 
gas are due to strong two-phase instabilities. See text for more 
discussion. Continued on Fig. |4] 
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Figure 4. Continuation of Fig. [3] 



Figure 6. Continuation of Fig. [5] The cold gas disk and the hot 
Mach disk are clearly seen at 2; = and z = —0.02, respectively. 
The hot 'plume' visible between z = 0.01 and 0.06 is actually a hot 
low-density sheath supported against radial collapse by rotation. 
Continued on Fig. [6] 











Figure 5. Edge-on view of the combined hot- wind-cold-disk gas 
evolution. Each panel describes the temperature in degrees K. 



Figure 7. (continuation of Fig. ^6^ 
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(a) isodensity surface 




(b) streamlines 

Figure 8. Density and streamline rendering of the gas flow at the 
same point shown in the last frame of Figs.[2}{7l Panel (a) show a 
level-set surface corresponding to gas density of 0.003 particles/cc. 
This surface is approximately cylindrically symmetric, but has 
been cut by a vertical plane passing through the axis of symmetric 
to reveal its structure: a hollow sheath with the shape of a hornet's 
nest beginning with the bow shock at the bottom which a large 
narrow downward dimple corresponding to the inner sheath. The 
surface is color coded by temperature, showing that the inner 
sheath is nearly 10^ K as it leaves the disk. Panel (b) shows 
the stream lines color coded by the temperature scale. The inner 
sheath in bounded by spiraling stream lines whose rotation comes 
from the angular momentum of the cold disk gas. 



300 Myr after the initial impact. The gas continues to flow 
around and ablate the cold disk. The post-encounter wind 
mixed with the ablated gas forms a hot sheath. Owing to the 
angular momentum of disk gas ablated from approximately 
3 disk scale lengths, the sheath has positive angular mo- 
mentum with respect to the cylindrical axis. The unshocked 
wind, flowing around the galaxy, pressurizes the sheath but 
the inward motion toward the symmetry axis is halted by 



Figure 9. Streamlines showing gas flow for a hot wind with an 
impact angle of 45 degrees with a cold gas disk at 373 Myr after 
impact. The gas density of the flow in units of atoms/cc is shown 
on the translucent plane. The overall morphology of the flow is 
qualitatively similar to that for face-on impact. 



its angular momentum barrier. This results in a relatively 
empty hot central core (final panel of Fig. [6|). 

By 300 Myr, the hot post-shock gas at impact param- 
eters smaller than 2 scale lengths has begun to push itself 
through the cold disk. This has three consequences: 1) the 
additional pressure provided by the hot gas drives a strong 
two-phase cooling instability; 2) as the hot gas is pushed 
through the hot holes created by the instability, it receives 
some of the angular momentum of the cold gas. This fiow 
is now partly rotationally supported against inward radial 
motion, creating a second inner sheath. Inside this sheath, 
there is only a hot tenuous gas with a much lower gas den- 
sity than that of the initial wind (by approximately an order 
of magnitude); 3) the shear created by the vertical oozing of 
the hot gas through the now clumpy cold-gas layer excites a 
Kelvin-Helmholtz instability in the disk plane (this will be 
described in detail below). 

The fully developed wind-cold disk interaction is de- 
scribed in Figures |4] and [71 The final panel of each of 
these figures describes the fiow at the end of the simulation 
(t ~ 500 Myr) . The hollow inner sheath begins to propagate 
upward in a spiral owing to the net angular momentum ob- 
tained from the disk. The inner sheath if very hot (T ^ 10^ 
K) but with low density (< 10^/cc). Figure [8] further illus- 
trates the state of the gas at 500 Myr. The upstream gas 
is gravitationally focused toward the galaxy giving the iso- 
density contours a hornet's nest shape. Panel (a) shows a 
rendered surface at a constant density of 0.002 atoms/cc 
color-coded by gas temperature in degrees Kelvin. At the 
bottom of the surface sits the Mach disk at a temperature 
of 8 X 10^ K. The projected rectangular feature between 
z = 0.1 and 0.3, x = —0.05 and 0.05 is the inner boundary 
of the sheath; you are seeing a cut through a low-density 
tube. Panel (b) depicts stream lines integrated through the 
gridded volume used to compute the ensemble field quan- 
tities, also coded by gas temperature. The location of the 
shock front is clearly seen by the sharp jump in the temper- 
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ature of the ambient wind gas from 10^ K to T > 4 x 10^ 
K. After impact, the stream hnes flare outward again and 
the gas cools a small amount. Approaching the inner galaxy, 
the hot gas attempts to push through the disk. In doing so, 
the hot gas gains some angular momentum and the cool gas 
loses some angular momentum. In the end, some net angu- 
lar momentum is lost to the hot gas and the resulting spin 
can be seen in the streamlines. A small fraction of the low- 
density gas squeezes through the multiphase disk at small 
radii and small specific angular momentum to form a funnel 
of hot gas with T > 3 x 10^ K. 

This same numerical experiment has been repeated with 
a wind direction of 45 degrees with very similar results. Fig- 
ure [9] describes the resulting flow with temperature-coded 
stream lines, as in Figure [8] with a density field slice along 
the symmetry axis. The stream lines and density slice re- 
veal the same inner and outer sheath morphology as Figure 
[9l The cold disk response and the angular momentum trans- 
port are also quite similar to the face-on case, to be described 
below, so we will not consider the details of the oblique- wind 
simulation any further. 

Typi cal of the galaxies wi th strong ICM interactions de- 
scribed in lChungetHI (|2009h . NGC 4522 shows little HI in 
the outer disk. Even more remarkable is the appearance of 
some gas above the inner disk, consistent with our predic- 
tion that the hot gas will force itself through the hot phase 
of the two-phase unstable disk during the later phases of 
the interaction. The extent of the neutral gas provides a 
further interesting diagnostic for the ICM-ISM interaction. 
Our simulation predicts that ram-pressure stripped hydro- 
gen will not be observable in its neutral phase (cf. Figs. 
[SHZl) owing to the heat provided by the shock and carried 
by the turbulent flow which rapidly heats the initially cold 
gas in the outer galaxy. There are some caveats: the ICM 
density used here is larger by an order of magnitude and 
the wind velocity used here is smaller by nearly an order of 
magnitude than those values appropriate for Virgo. Given 
the complexity of this interaction, a direct prediction should 
be made with the appropriate values. In addition, because 
this simulation does not include gravitation on small scales; 
it is possible that this simulation misses small scale cold- 
blobs forming from cooling instabilities. However, the per- 
vasive shock heating suggests to me that HI tails detected 
in galaxie s that also show evidence of r am-pressure interac- 
tions (e.g. IVollmer k, HuchtmeierlbOOTI ) have lost their neu- 
tral gas through tidal stripping. Nonetheless, the morphol- 
ogy of gas in NGC 4522 is qualitatively similar to that seen 
in these simulations and certainly ap pears to be an o ngoin g 
wind-disk interaction as indicated bv lVollmer et al.l (|2008l V 
It is conceivable that both mass-loss mechanisms are occur- 
ring simultaneously in NGC 4330 (Abramson et al. 2011). 
The time scale for survival of tidally-lost neutral gas in the 
hot turbulent wake of the ICM-ISM interaction may be dis- 
tinctly morphology dependent and will require further work. 
Future work using the full DSMC collisional implementation 
will yield detailed radiative diagnostics (see 3S1) which will 
help disentangle the mechanisms. 

Figure [10] shows the azimuthally averaged distribution 
the disk's gas mass inside of a given radius (in kpc) and as a 
function of time (in Myr). We see that approximately 40% 
of the initial gas mass lost, and this is mostly from the outer 
disk and in the first 300 Myr. As seen in Figures [2HZ1 the 



Gas fraction 




Figure 10. The cumulative fraction of gas in the disk with radius 
(in kpc) and time (in millions of years). 



Cumulative angular momentum 




Figure 11. Similar to Fig. 1101 this figure shows the evolution 
of the cumulative angular momentum fraction of gas in the disk 
with time. The contours on the base plane show fixed fractions 
for interpretive guidance. 



initially cold gas in the outer disk is ablated and heated. 
Although the gas continues to be ablated at a low level, the 
loss rate is diminishes rapidly after 300 Myr owing to both 
the increased density of the cold gas after pressurization and 
by the decrease in the scale length of the cold gas after its 
sizable loss of angular momentum. 

Figure [11] shows the evolution of the angular momen- 
tum fraction with time. This further emphasizes that the 
disk gas is peeled off from the outside in and demonstrates 
that the overall gas content of the inner galaxy is largely un- 
affected. More than 60% of the original angular momentum 
content of the cold gas disk is lost to the wind, the dom- 
inant fraction owing to the mass loss. The wiggles in the 
constant angular momentum fraction contours are due to 
planar Kelvin-Helmholtz instabilities which we will describe 
next. 

Figure [12] shows a sequence of density cuts through the 
disk plane. The density scale on the left of each panel de- 
scribes the color code for density in the range n G [1.5 x 
10~^, 1.5 X 10^] atoms/cc. The hot gas pushing through the 
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Figure 12. Gas density as a function of time (in system units). System density units are 1.5 x 10 ^ atoms/cm^ and one scale length is 
0.01 system length units on the x and y axes. The time at the top of each frame is shown in millions of years (Myr). 



cold disk from below picks up a small amount of angular 
momentum. This, in turn, creates a strong shear layer in 
the disk plane. This shear gives rise to a rampant Kelvin- 
Helmholtz (KH) instability. The KH waves develop and 
break, giving the familiar saw tooth pattern. The toothed 
ring then fragments: the inner parcels continue to move in- 
ward and the outer parcels move outward consistent with 
their specific angular momentum values. Gas from the de- 
stroyed ring coalesces into new rings and the process repeats. 
This limit-cycle behavior is seen as saw-tooth-like wiggles in 
the angular momentum contours in Figure [TT] At approx- 
imately 200 Myr into the wind- impact event, a stellar bar 
forms. The influence of the the bar on the flow pattern is 
noticeable at T ^ 350 Myr. 

The extent of the multiphase medium is shown in Fig- 
ure [13] for the disk plane. Panel (a), and the halo. Panel 



(b). This density-temperatur e plot shows the classically ex- 
pected two-phase instability (|Field et al.lll969l ) . The hot gas 
compresses the initially cold in the inner disk and shortened 
cooling time keeps it cold. The initially cold gas in the outer 
disk is warmed by the turbulent interaction with hot wind 
and this gas create the outer sheath. This results in the con- 
tinuous distribution between the two phases in the figure. 
Note that this transitional gas only exists in the vicinity of 
the disk. Panel (a); the sheath itself has T ^ 10^ K, Panel 
(b). Prior to the wind interaction, the T-p phase plane shows 
the hot tenuous wind and cold disk gas in two distinct and 
separated modes. 

Figures [15] and [16] show the density and temperature 
for the gas at 100 Myr from the sequence in Figure [12] It 
is a particular moment in the non-linear evolution of the 
KH instability excited by shear created by the low angu- 
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Figure 13. Density-temperature phase plane at 200 Myr for the 
initially cold disk gas (top), confined to the region of the stellar disk, 
and the initially hot wind gas (bottom), filling the region initially 
below but later above and below the stellar disk. The logarithmic 
color scale denotes relative density on the phase place. The initially 
cold gas develops a distinct hot tenuous phase through interaction 
with the hot wind (top). Similarly, the hot wind is both heated and 
cooled through its interaction with the cold disk gas (bottom). 
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Figure 14. As in Fig. 1131 but for the simulation including hot 
coronal gas initially. 



lar momentum gas forcing its way through the disk. This 
particular snapshot shows a combination of coalesced rings 
connected by the spokes of the remnant KH instability. A 
detailed check shows that the hot and cold phases are in ap- 
proximate pressure equilibrium. Moreover, the classic sat- 
urated KH instabilities are clearly seen on multiple scales 
(i.e. examine the saw-tooth shaped interface between the 
high and low density gas in these figures). Unlike previous 
figures that render ensemble- averaged field quantities, this 



figure shows individual gas simulation particles color coded 
by their ensemble field values. 

Finally, we may ask, are such things seen in Na- 
ture? I would argue that they are. For example, compare 
NASA/ESO image of the Cartwheel galax^fj with the simu- 
lation profiles at 200 Myr fFig. [T7|) . Commonly, these rings 
and spokes are interpreted as a splash created by an in- 
truder. A close look Cartwheel image galaxy reveals the 
same spokes and KH scalloping seen in Figure [17] suggest- 
ing that the shear between hot wind gas and the cold disk 
gas is shaping the in-plane features! ROSAT and Chandra 
observations of the Cartwheel system reveal a stream of soft- 
x-ray emitting gas co nnecting the Cartwheel to its apparent 
companion intruders (|Wolter fc Trinchierill2003l l with an in- 
ferred column density of Nh = 1.3 x 10^^ cm~^. Assuming 
an outer ring diameter of 60 kpc, one estimates a hot gas 
space density of roughly uh = 5 x 10 ~^ cm~^, comparable 
to t he wind densities in the simulations here. 

Higdon (l996) estimated the expansion rate of the ring 
to be 53 =b 9 km/s based neutral H observations, consistent 
with the passage of a n intruder approximately 300 Myr ago. 
lAmram et al.l (| 19981 ) estimated a mean expansion velocity of 
13-30 di 10 km/s with a superimposed sinusoidal distortion 
of amplitude 20 ± 5 km/s from Ha observations. However, 
the expansion rate of the ring is similar to that found in 
the non-linear KH instability due to the wind. Altogether, 
it seems plausible that the observed spokes are not, as some 
have suggested, simply the original spiral arms of the target 
galaxy reappearing but the same KH instability seen in the 
simulations described here. 



3.5 An edge-on Avind 

The affect of the edge-on wind differs qualitatively and quan- 
titatively from that of the face-on wind. First, as might be 
expected, the shock forms in a rim around the outer cold 
disk, not like a Mach disk. The post-shock gas mixes with 
and ablates the cold gas. The angular momentum of the ini- 
tially cold gas gives the flow some spin relative to the disk. 
The mixture of incoming wind and ablated cold gas moves 
around the disk in the direction of disk rotation. Approx- 
imately 70% of the cold disk's initial angular momentum 
content is lost to ablation in this way. However, only ap- 
proximately 10% of the cold disk mass is lost by 500 Myr 
after the encounter; this is nearly all of the gas outside of a 
radius 0.04 (12 Kpc in Milky Way units). Unlike the result- 
ing disk in face-on encounter, there is little shrinkage of the 
inner disk. 

At the same time, the hot wind flows under and over 
the disk and subsequently this gas is gravitationally focused 
on to the face of the disk. This, in turn, pressurizes the cold 
disk and drives a strong multiphase instability. The pres- 
sure from the gravitationally focused wind is higher than 
from the Mach-disk-protected face-on wind leading to higher 
density in the cold gas phase. The phase-plane distribu- 
tion is similar to that for the face-on wind but very little 
of the gas has transitional densities and temperatures; the 
two phase remain well separated. In summary, although the 
overall outcome is similar to that of the face-on interaction, 



^ e.g. http : //chandra . harvard . edu/photo/2006/cartwheel/ | 
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Figure 15. Cold disk from the simulation from Fig. [12] at 100 Myr. The left-hand panel shows gas density and the right-hand panel 
shows temperature. 
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Figure 16. Cold disk at 100 Myr: close up. Panels and scales as in Fig. 1151 



many features of the dynamical pathway to the final state 
are different. 



3.6 The addition of coronal halo gas 

To get some insight for the overall change in the wind-cold- 
disk interaction caused by a hot coronal halo, we repeat 
the simulation described in detail in ^3.41 but adding a hot 
isothermal coronal component in hydrostatic equilibrium in 
the galaxy with an initial temperature of 10^ K and a density 
of 10~^ atom/cc at a radius of 30 kpc. This value was cho- 
sen to match typical hot halo densities estimat ed at large 
galactocentric radii. However, as described by iFang et alj 
(|2013|), this value may overestimate the inner galaxy den- 



sity and underestimate the total coronal gas mass owing to 
the assumption of a cuspy, NFW profile. 

The main results are as follows. The pressure from the 
hot- wind impact is more readily communicated through the 
coronal halo that is already in place downstream of the 
shock. This increases the pressure on the inner disk and 
results in a higher amplitude multiphase instability. This is 
seen through comparison of the density-temperature phase 
plane at the same interaction time with and without a 
corona (Figs. [13] and [14). Recall that the upper panel shows 
the gas in the disk and the lower panel shows the gas in the 
halo. Overall one sees more high density cold gas in the sim- 
ulation with a corona. In addition, one sees that the there is 
smoother transition between the cold-dense and hot-tenuous 
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Figure 17. Cold disk at 200 Myr. Panels and scales as in Fig. [TS] Note that the hot gas sits in the void opened up by shear in the 
cold gas. The saturation of the KH instability gives rise to rings and spokes, similar to those seen in ring galaxies. Moreover a close 
examination of the ring in the Cartwheel reveals the same sort of scalloping seen in the simulation shown here. In both panels, the scalar 
is logarithmically mapped to a color sequence that transitions from blue to red to yellow to white. On the left, the density of the hot gas 
(red) is 10~^ cm~^ and the density of the cold gas is 10^ cm~^. On the right, the host gas (yellow- white) has T k, 10^ K and the cold 
gas (blue) has T ?^ 5000 K. 



phases in the flow without a corona; the corona suppresses 
the production of the intermediate-state gas by providing an 
additional heat source. 

Similarly, the higher density of the hot gas at smaller 
galactocentric radii increases the angular momentum trans- 
fer from the cold disk to the post-interaction outflow by 
increasing the collisions of the hot, low angular momentum 
gas and the cool, high angular momentum gas. This is seen 
in the comparison of the cumulative angular distributions 
shown in Figures [TT] and 1211 The total angular momentum 
drops nearly linearly to 10% of its initial value at 450 Myr 
after the initial impact. However, at the same time, the re- 
moval of this angular momentum causes the cold gas disk to 
shrink. This decreases the overall mass loss. In other words, 
the net result of the coronal gas is to protect the initially 
cold gas mass against stripping through ablation (compare 
Figs. [TOl and [20]) . The disk loses only approximately 20% of 
its total mass overall, but by 450 Myr, all of this gas is inside 
of 8 kpc. 

The overall flow pattern is the same with the hot corona 
(Fig. [19} and without (Fig. [H lower panel). The extent of 
both the inner and outer density sheaths are smaller with the 
corona, consistent with larger angular momentum loss and 
disk shrinkage. Moreover, the simulation clearly shows that 
the density of the cool gas in the inner disk is much higher 
than without the coronal gas. The effect of the coronal gas on 
the inner galaxy may itself place a limit on the density of the 
inner coronal gas and thereby also place a limit on the shape 
of the galactic dark matter potential. This scenario is ripe 
for further investigation. Use of the non-LTE DSMC code, 
under development, will allow prediction of line- strengths for 
observed density and temperature diagnostic features (e.g. 
Mg II absorption). 



4 DISCUSSION 

The effect of ram pressure on galaxy morphology and trans- 
formation has been studied in several contexts over the 
last decade. A pa rticularly nice example is provided by 
ICrowl et al.l (|2005[ ) who inferred the signature of ram pres- 
sure and ablation by combining optical, HI, and radio con- 
tinuum observatio ns for NGC 4402 (cf. Figs. 2 and 3 from 
ICrowl et al.ll2005[ ). They argue that the absence of redden- 
ing at the disk's edge implies truncation and that the po- 
sition angle of linear dust filaments coincides with radio 
continuum, su^^estin ^ the ICM wind orientation. Similarly, 
IVollmer et al.l (|2008|) show that observed regions of extra- 
planar gas with H-a flux have larger velocity than those 
with CO flux in Virgo galaxy NGC 4438, suggesting that 
the ionized gas is excited by the ram-pressure mechanism. 
All in all, there is little doubt that ICM winds influence a 
galaxy's ISM. 

These and related observations have motivated a variety 
of computational hydrodynamic efforts and provide oppor- 
tunities for comparison. The morphology of the flows in our 
simulations have the same morphological features as t hose 
performed with hydrodynamic codes (e.g. iBekkil |2009|, see 
his Figs. 2 and 7, using GRAPE-SPH). We expect DSMC to 
have lower resolution for the particle numbers used here than 
many of the hydrodynamic simulations discussed above, and 
therefore the good correspondence is encouraging. Also as 
expected, the interfaces between gas components and multi- 
ple phases are distinct, demonstrating a virtue of the kinetic 
approach. 

However, since most authors had different sets of inter- 
ests, goals, and physical choices, direct intercomparisons are 
often approximate. In addition to different initial conditions, 
many used different equation of state combinations with 
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Figure 18. Streamlines showing gas flow for an edge-on impact 
of a hot wind with a cold gas disk at 400 Myr after impact. The 
wind flows from left to right in this flgure and is color coded by 
gas temperature (K). The gas density of the cold disk layer with 
density of 1 atom/cc is shown as the white surface. The up-stream 
gas is gravitationally focused toward the galaxy. The gas shocks 
on impact and the hot flow ablates the cold gas as in the face-on 
impact. Here, the ablated gas is pushed over or under the disk; 
through mixing, the mixture of the wind gas with disk gas gives 
the flow some overall rotation. Some of the hot gas is refocused 
down to the disk plane, heating the disk from above. 




Figure 19. Similar to Fig. [S] streamlines showing gas flow for an 
edge-on impact of a hot wind with a cold gas disk at 250 Myr after 
impact. The flow is qualitatively similar to that without a corona. 
The main difference is the higher temperature of the downstream 
flow at small cylindrical radius. 
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Figure 20. The cumulative fraction of gas in the disk with radius 
(in kpc) and time (in millions of years) for simulation with a 
coronal halo (cf. Fig. I10|) . 



Cumulative angular momentum 




Figure 21. Evolution of the cumulative angular momentum frac- 
tion of gas in the disk with time for simulation with a coronal halo 
(cf. Fig.[ll]). 



and without atomic an d plasma cooling. For some exam- 
ples, iQuilis et alj (|2000|) simulated the interaction between 
a spiral galaxy and the ICM of a rich cluster at a relative 
speed of 2000 km/s. They show that the interaction removes 
most of the gas in 100 Myr and argue that cooling is not im- 
portant in this case. Our wind speeds are nearly an order of 
magnitude smaller than this, yet encouragingly, the morpho- 
logical features of the ram-pressure response are very similar 
to those found here (see their Fig. 1). In contradiction, we 
conclude that the wind-disk interaction cannot remove all of 
the initially cold gas. However, we have not tried such a fast 
encounter; given the greatly increased energy and momen- 
tum content in this fast wind, it see ms plausib le that our 
result would be sim ilar to theirs. Marcolini et al.l (|2003) and 
iMaver et al.l (|2006h performed hydrodynamic simulations of 
ram pressure stripping for dwarf galaxies. The former inves- 
tigation considers winds typical of poor galaxy groups and 
considers the role of mass, velocity and wind direction. They 
find that dwarf galaxies may be stripped in a few hundreds 
of Myr (see their Fig. 4 for an illustration of the interaction 
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morphology). Mayer et al. (2006) pattern their study on a 
Milky- Way-hke system and consider the com bined roles of 
ram p ressure and orbitally induced tidal forces. iBoselli et al.l 
(|2009!) attempted galaxy evolution modeling with ram pres- 
sure and starvation. Their work suggests that anemic spirals 
with truncated gas distributions and young populations have 
produced by recent ram pressure event and cannot be pro- 
duced by starvation alone. Similarly, they report that low 
luminosity dwarf disk galaxies can be transformed in dE 
galaxies with ease. Most of these are not directly compara- 
ble to the simulations described in ^3.4l but they still exhibit 
many of the same morphological features qualitatively. 

Many researchers compare their stripping results to the 
predictions of the Gunn-Gott criterion, although the gen- 
eral consensus is that this is only a crude approximation. In 
particular, as in this paper, a number of hydrodynamic simu- 
lations indicate an early phase of rapid stripping in the outer 
disk followed by a period of slower, viscous, str ipping (e.g. 
iMarcolini et al.ll2003l : [Roe diger fc Henslerl 12005V However, 
we demonstrate here that the initial mode of rapid stripping 
is ablation, not impulsive mo mentum transf er. Based on an 
extensive parameter survey, iRoediger fc Hensleri (|2005i ) de- 
velop a revised criterion based on pressure at the disk mid 
plane. However, it is clear that the mass loss in the end de- 
pends on the widely diffe ring ICM densities, velocities and 
galaxy mass (|Bekkill2o69l V 

In this study, our goal was the interaction physics itself 
rather than the implications for a specific astronomical sce- 
nario. In summary, although momentum is transferred from 
the wind to the disk ISM, the process bears little resem- 
blance to its moniker, ram pressure, but rather seems to be 
an elaborate interplay between shocks, shear fronts, kinetic 
ablation (i.e. evaporation) and angular momentum transfer. 
Although it is natural to seek a scaling relation such as the 
Gutt-Gott criterion to summarize the overall outcome, the 
results of this paper suggest that an impulsive ram pressure 
approximation will be crude at best. The large variation of 
successes and failures of the criterion reported in the litera- 
ture is not a surprise in this light. 

Motivated by he discovery of a large (M ^ 3x 10^ Mq) 
HI cloud in the Virgo cl uster that appears to origina te 
from the galaxy NGC 4388 (jOosterlo o k van Gorkom' 2005) , 
iRoediger et al.l (J2006) used hydrodynamic simulations to 
confirm that ram pressure is sufficient to strip this galaxy's 
gas mass. They report a distinctive flow pattern that they 
ascribe to von Karman oscillation, typically produced by 
the unst eady flow of a fluid past a blunt body, and suggest 
that the lOosterloo fc van GorkomI (|2005l ) discovery is evi- 
dence of this. The authors do remark that the observed HI 
tail morphology does not match the simulation and suggest 
that this discrepancy is due to missing physics and ICM 
interaction. We agree that this is likely. In particular, the 
results of t his paper challenge th e assumption that an adi- 
abatic gas ([Roediger et al.l 120061 ) will correctly capture the 
nature of the mesoscopic interactions important at interfaces 
and that the hydrogen lost fro m the stripped galaxy will be 
neutral. As in the present work jRoediger &: BriiggenI (|2006|) 
found that wind orientation angle does strongly affect mass 
loss unless the wind direct ion is nearly edge on. 

iTonnesen et al.l ([200?) studied the problem from the 
ACDM galaxy formation and evolution viewpoint by track- 
ing galaxies in clusters over time in a cosmological volume 



to ascertain the distribution of evolutionary drivers. They 
conclude that gas stripping is the most common mass loss 
mechanism and remains important out to the virial radius of 
the cluster. This general finding is consistent with our sim- 
ulations. However, their aim was an ensemble survey using 
ACDM initial conditions. Given their limitations on resolu- 
tion, a direct comparison with the details of the dominant 
physical mec hanisms in thei r simulations is not possible. 

Finally, Kapf erer et al.l ([20091 ) explicitly investigated 
the dependence of star formation and morphology on the 
strength of ram pressure using combined N-body and hy- 
drodynamic simulations with recipes for star formation and 
stellar feedback. They find that star formation is enhanced 
by more than a magnitude in the simulation with a high 
ram pressure relative to the same system evolving in isola- 
tion. They report that newly formed stars in the disturbed 
gas stars fall back to the old stellar disc, building up a the 
bulge. They find also that gas stripping saturates quickly as 
in the simulations reported here. In addition, they describe 
a complex velocity pattern due to the rotation and spiral 
arms. This is qualitatively identical to the KH instabilities 
described in ^3.41 

All in all, the simulations reported in this paper describe 
similar features and mass loss fractions to those from hydro- 
dynamic simulations when direct comparisons are possible. 
Moreover, our simulations underscore the importance of the 
interface in predicting the evolutionary morphology of ram- 
pressure induced flows and response. This suggests that the 
kinetic approach may become an essential tool for clarifying 
the dynamics of multiscale gas interactions for astrophysical 
flows. 



5 SUMMARY 

The ram pressure mechanism creates flows in multiply dis- 
tinct gas phases with dynamically important interfaces, 
such as shock fronts and phase boundaries. Resolving the 
physics at these interfaces requires subtlety when using a 
hydrodynamic method but is naturally done by direct solu- 
tion of the collisional Boltzmann equation using a kinetic- 
theory method. To this end, we explores the use of the hy- 
brid kinetic theory-n-body code from Paper 1 to provide 
a check on our physical understanding of the ram-pressure 
scenario. The kinetic theory algorithm. Direct Simulation 
Monte Carlo (DSMC), is computationally costly, so we have 
not extensively surveyed the range of parameters relevant 
to all ICM-ISM interactions but rather focused on the dy- 
namical mechanism for a single case, more typical of a 
group ICM density and vel ocity for a Milky- Way like galaxy 
([Freeland &: Wilcotsl [20T1I ) . Overah, the findings from the 
simulations here are consistent with many of the previously 
published results. The simulations include a gravitationally 
self- consistent stellar disk, spheroid, and dark matter halo 
with an initially cold gas layer. The wind temperature is 10^ 
K with a uniform bulk velocity set to twice the virial velocity 
of the galaxy initially. The choice of a group rather than rich 
cluster environment is partly motivated by computational 
expense. The key physical features of the gas response are 
as follows: 

(i) The incoming wind shocks as it approaches the disk 
and a compression wave begins to propagate upstream. 
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Meanwhile, the post-shock gas is heated, compressed, and 
this results in a lower bulk velocity in the downstream gas 
by approximately a factor of four. The ballistic crossing time 
of the hot gas is now close to the disk crossing time. The 
velocity of the hot wind gas decreases further as it pushes 
its way through the multiphase ISM. 

(ii) The cold gas is not launch ed from disk en ma sse as 
in the impulsive approximation (|Gunn &; Gottlll972l V but 
rather, the hot post-shock gas ablates and entrains the ini- 
tially cool disk gas. The stripped gas ends up with a tem- 
perature close to that of the entraining wind. 

(iii) The hot wind, either pushing through the disk for a 
face-on and oblique flow or impacting the disk from below 
and above for an edge-on flow, pressurizes the cold disk and 
promotes a two-phase instability. The rapid consumption 
of gas by star formation and removal of the cold gas supply 
may well result in morphological transformations as pre vious 
authors have suggested. Indeed, iKapferer et al.l (|2009|) and 
iTonnesen fc Brvan (2010) find that gravothermal instability 
may promote cooling and star formation; our simulations 
neither include the molecular cooling or gravity on molecular 
cloud scales needed to investigate this possibility. 

(iv) In the face-on or oblique wind-disk encounters, the 
post-shock ICM flowing through the now porous ISM has 
low specific angular momentum relative to the rotationally 
supported ISM. This results in a strong shear between the 
two phases with two important consequences: 1) the shear 
excites a strong Kelvin-Helmholtz instability that results in 
"ring and spoke" morphology; and 2) the hot wind carries 
off some the gas disk's angular momentum. This results in 
an inward advection of cold gas with the possible formation 
of an inner ring and a net rotation of the wind escaping 
the disk. The addition of angular momentum to the coro- 
nal plasma with the rotational orientation of disk may have 
interesting implications for galactic dynamos. 

(v) Since the escaping gas is primarily lost from the outer 
disk, the outflow appears as a bell-shaped sheath for the 
face-on and oblique flows. The angular momentum gained 
through the interaction appears to help collimate the sheath 
by providing a centrifugal barrier at small radii. For an edge- 
on encounter, there is no well-defined inner sheath, although 
the cold gas disk casts a geometric shadow on the down 
stream flow creating a density dip along the flow axis. 

In summary, we find that some initially cold gas does 
leave the galaxy as do most studies, but our simulations 
suggest that the process is complex. It is, therefore, unlikely 
be shoehorned into a simple scaling law owing to the sub- 
tle dependence on initial conditions. The work here further 
suggests that the gas leaving the galaxy due to the wind 
interaction will be hot (i.e. not detectable as neutral hy- 
drogen) . We would need higher resolution and huge particle 
numbers to see the possible entrainment of small cold clumps 
or to resolve gravothermally unstable clumps. Furthermore, 
these simulations show that radiative cooling is critical to 
the flow pattern and cannot be ignored. Without cooling, 
most of the cold gas would heat and subsequently evapo- 
rate whereas with the cooling, we find strong multiphase 
instabilities that protect the gas against further large-scale 
mass loss. More generally, the phase and shock interfaces are 
likely to control the overall evolution of the gas and produce 
structure at all scales! Similar to collisionless dynamics, one 



can not gloss over interface regions even if their total volume 
is small. 

Although DSMC was developed for high-Knudsen num- 
ber (i.e. rarefied) flows, it may be used for near continuum 
flows as well. DSMC is more computationally expensive than 
traditional hydrodynamics for near-continuum flows, but 
has the advantage of unconditional stability, self- consistent 
momentum transport and no difficulties with formal dis- 
continuitiecl- In particular, an ideally tuned DSMC simu- 
lation simultaneously maintains approximately 10 particles 
per simulation volume whose length scale is approximately 1 
mean-free path. For denser flows, this condition is computa- 
tionally infeasible. The implementation described in Paper 1 
has a collisional limiter that transitions to the near equilib- 
rium solution of the energy and momentum flux equations 
when the collision rate becomes very high. However, even 
in this limit, momentum and energy remain precisely con- 
served and numerics of the flow remains stable. This makes 
DSMC ideally suited for investigating flows with momentum 
transfer and heating or cooling by conduction or radiation 
at interfaces. 

Since this project arose as a code test, we have chosen to 
use an LTE approximation for gas cooling; this approach is 
also known as the quasi static state (QSS) approximation. 
QSS describes that collisional and radiative excitations as 
rates based on the local fleld quantities. In addition, we do 
not include electrons as a separate species for these simula- 
tions, but include their influence in the QSS limit. There- 
fore, electron conduction is not self-consistently simulated 
although hard-sphere conduction of the atoms is. On the 
other hand, a precise study of the excitations and electron 
interactions at the interfaces is necessary for estimating their 
physical state and dynamical influence on the flow. More- 
over, accurate observational predictions depend on the state 
of the gas at the interface. It is precisely such investigation 
of interface dynamics and energetics on small and intermedi- 
ate scales that motivated the combined DSMC-n-body code 
for transitional flows in the flrst place. 

Therefore to improve the self consistency of these sim- 
ulations, we are actively developing a version of our hy- 
brid code that inc l udes the full CHIA NTJj atomic database 
(|Dere et al.l [1997; Lan di et al.l 1201 ih for collisional cross 
sections and recombination (bound-bound and bound-free) 
and standard plasma transitions (free-free) for the ionized 
regime. This new implementation will allow us to straight- 
forwardly make explicit predictions for both continuum and 
line-excitation strengths. However, it does not yet include 
the electrons as a separate kinetic species, but rather as- 
sumes that they follow the ions. This restriction will be re- 
laxed in a future version that will allow a more realistic 
treatment of electron conduction and the dynamical influ- 
ence of flxed magnetic flelds. To be sure, simulation of these 
processes based on kinetic theory are very expensive but, in 



^ In this sense, one might say that DSMC is 100% shock cap- 
turing, although more precisely DSMC has no need to capture a 
discontinuity because the physics of a shock front is resolved by 
the algorithm at the scale of the effective mean free path. 
^ CHIANTI is a collaborative project involving George Mason 
University, the University of Michigan (USA) and the University 
of Cambridge (UK). 
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principle, will yield detailed observational predictions and 
accurate internal energetics and momentum transfer. 
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